Inhomogeneous LOFF phase revisited for surface superconductivity 
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We consider 2D surface superconductivity in liigli magnetic fields parallel to the surface. We 
demonstrate that the spin-orbit interaction at the surface changes the properties of the inhomoge- 
neous superconducting Larkin-Ovchinnikov-Fulde-Ferrell state that develops above fields given by 
the paramagnetic criterion. Strong spin-orbit interaction significantly broadens the range of exis- 
tence of the LOFF phase, which takes the form of periodic superconducting stripes running along 
the field direction on the surface, leading to the anisotropy of its properties. In connection with 
experiments by J.H. Schon et al. [Nature 914, 434 (2001)] on superconductivity of electrically 
doped films of the cuprate material CaCu02 , we also discuss this problem for the d-wave pairing to 
indicate the possibility of a re-orientation transition as the magnetic field direction is rotated in the 
plane parallel to the surface. Our results provide a tool for studying surface superconductivity as a 
function of doping. 
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There has been renewed experimental and theoreti- 
cal interest in the properties of metallic states localized 
at a surface, and sm'face superconductivity (SSC). Sur- 
face states (a.k.a. Tamm's levels) are well-known from 
the physics of semiconductors. Numerous ARPESjl| |^ 
data and STM studies of Friedel oscillations]^ have now 
proven the existence of 2D metallic bands even at sur- 
faces of metals. The bands are well-separated from the 
bulk and possess clear-cut 2D Fermi surfaces. For heavy 
enough elements such Fermi surfaces are split further by 
strong spin-orbit (SO) interactions. For example, SO en- 
ergy for electrons at the Fermi level for Au is estimated 
at 0.1 eV 1^, while for Li-doped surfaces of Mo and W its 
value increases up to 0.13 eV and 0.5 eV, correspondingly 
§. 

Islands of a surface superconducting phase were 
also observed for the surface-doped tungsten bronzes, 
WOsiNa, at = 91.5i^@ . The ARPES resuhs men- 
tioned above suggest that SSC may actually be a rather 
widespread phenomenon. 

It is noteworthy that the bulk WO3 is an insulator 
at low doping[H. Recently a remarkable breakthrough 
has been accomplished by doping insulator surfaces elec- 
trically in the so-called Field-Effect Transistor (FET) 
geometry Q. SSC with high Tc was induced by both elec- 
tronic and hole doping of films of the prototype cuprate 
material, CaCuO2[0|. 

All this makes us believe that the search for SSC 
emerges as a new and important development in stud- 
ies of the properties of surfaces, especially their metallic 
properties. As for the SSC itself, its mechanisms are un- 
known, and may have nothing to do with the ones in the 
bulk. Of a special challenge is the possibility of super- 
conductivity at surfaces of ordinary metals such as Cu 
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Q or lithimn-doped Mo and WQ, and the influence of 
substrate on superconductivity in thin films. Thus, it 
is necessary to return to a more careful investigation of 
the low temperature properties of surfaces with absorbed 
atoms and the role of adsorbed atoms as dopants of car- 
riers into the metallic surface zone (e.g., see Ref.Q) like 
it was discovered for WOatNa 

The major challenge to such a program lies in discern- 
ing SSC and studying its properties. Even for insulators 
where doping by FET, ideally, provides an effective con- 
trol on surface properties, experimental tools to probe 
SSC are limited in numbers. The Meissner effect due 
to surface superconducting islands would probably never 
produce a bulk screening. Thermodynamical probes are 
also difficult because of the smallness of contributions 
from surface layers of atomic thicknesses. So far, e.g., 
in Ref. surface superconductivity has only been de- 
tected by measuring resistivity dependence on tempera- 
ture in the perpendicular-to-plane magnetic fields. 

In what follows we focus on destroying SSC by mag- 
netic fields applied parallel to the surface. In high enough 
fields one should expect the appearance of a 2D version of 
inhomogeneous superconducting state known as Larkin- 
Ovchinnikov-Fulde-Ferrell phase (LOFF) It was also 
shown that two-dimensionality broadens the region of the 
LOFF state on the B-T phase diagram]^. Motivated by 
experimental findings mentioned above, we investigate 
the peculiarities introduced into this phenomenon by the 
SO effect or non-s-wave pairing. Experiments by Schon 
et al. not only have proved FET to be an effective dop- 
ing process. Different levels of doping result in differ- 
ent Tc-s. Our results for SSC in parallel magnetic fields 
are expressed in terms of this Tc, providing the tool for 
comparing theoretical predictions with experiments by 
controlling the doping level. Theoretically, there is no 
long range order in a 2D superconductor. However, cor- 
relations are destroyed on the exponentially large spatial 
scale, R ~ £,oexp{Ep /T), which would exceed the size of 
the film. 
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We employ below the weak-coupling BCS-like scheme 
by assuming that electrons interact via a weak short- 
ranged interaction, U{r,r'). Then Tc ^ ep, and only a 
narrow vicinity of the Fermi surface is involved. Thus, 
the interaction in the momentum representation can be 
taken in the form: 



(1) 



where p, p' lie on the Fermi surface; the angular de- 
pendence is expressed through a complete set of basis 
functions xi (p) (index / enumerates different representa- 
tions, as in expansions over the spherical functions in a 
3D isotropic model). Superconducting order parameter, 
the "gap", Aq,^(p) is defined by the equation: 



Ao/3(p) = \Ui\xi{p) 



where Fap{p']iujn) stands for the Fourier component of 
Gor'kov anomalous function: 



) = -(Tj*„(r,r)#Mr',r') ), (3) 



and Ui < is a constant in Eq.(|l|) corresponding to the 
selected pairing channel. When the field operators for 
electrons are re- written in momentum space, (r, r) — 
E„*a(p,T)e*P-, Eq.d) becomes: 



Aa/3(p) = |f/i|X/(p) 



0^Xiip'){^c.{p',r)^p{-p',r)). 



^ (4) 

The two operators inside brackets in Eq.(0) anti- 
commute. In the presence of the center of inversion (CI) 
the behavior of xi{p') at p' — > —p' alone determines the 
symmetry of Aq,^(p) (even (singlet) vs odd (triplet) par- 
ity pairing). For an s- or d-pairing the order parameter 
below Tr has the form: 



Aa/3(p;9) = A(q,T)(icrj^)c,/5Xi(p), 



(5) 



where the momentum q stands for the spatial dependence 
of the gap amplitude, A(q, T). 

Surface always breaks the CI symmetry due to the dif- 
ference between the " top" and the " bottom" . The direc- 
tion bulk-to-surface determines n, a unit vector normal 
to the surface. Qualitative changes in the surface elec- 
tronic spectrum come about from the well-known Rashba 
term@: 



hso = a{a- xp-n), 



(6) 



which specifies SO interactions at a surface. Eq.(|^) lifts 
the two-fold spin degeneracy for band electrons. The 
electron spectrum now consists of two branches with two 
Fermi surfaces: 

a 

e±{p) = vf{p - PF±); PF±^PFil± — ) (7) 

VF 



Tc 
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FIG. 1; Calculated Tc as a function of magnetic field and 
spin-orbit interaction app for a 2D surface superconductor. 
All values are in units of Tco, Tc for the 2D superconductor in 
the absence of magnetic field. Thus, B (ibB/Tco, hso 
app/Tco, and ^ Tc/Tco- 



Even though SO splitting, 2apF, may be on a scale of 
tenths of eV |l|, |, we assume that 2apF ^ e_F. Below 
we speak of a strong or weak SO meaning the relative 
values of 2apF and Tc. 

Some theory issues regarding superconductivity with- 
out CI due to the presence of the SO term Eq.(||) were 
first considered in Ref. and more recently in Ref. . 
The Gor'kov function (in the momentum representation), 

Fap{p,+0) = —{^a{p)^i3{—p)), that stands under inte- 
gral in Eq.(^) represents the wave function for Cooper 
pairs in the condensate. In the presence of CI symmetry 
the latter can be classified according to the parity: 



i^a/3(P;+0) = 



i{(Ty)ai}f{p); {f{p) even) 
i[{d{p) ■ &)ay]ap; (d(p) odd) 



(8) 



With CI broken by non-zero SO term Eq.(|g), the pairing 
wave function becomes a mixture of even and odd terms. 
It is important to realize that while this mixing changes 
physical properties of the SC phase, the gap order param- 
eter, Aa/iijj; q), preserves its singlet form Eq.(H). For in- 
stance, s-pairing indeed induces the non-zero triplet com- 
ponent Eq.(||), as shown in Refs.[^, However, the 
latter does not automatically generate a "triplet" gap, 
A^^(p;q). Indeed, re- writing the integration over p' in 
Eq.(^ as d'^p' dS'pd^, where ^ = vpip — Pf), we 
notice that the triplet _F-component is odd in particle- 
hole transformation, ^ — > — ^, and, hence, the integrals 
of the form Eq.(^ would only give small terms of or- 
der (apF/eF) <C 1. In other words, while SO interaction 
may significantly change spin structure of the normal and 
anomalous Green functions, the "gap" , Aaf3{p, q), Eq.(||) 
preserves its usual form Eq.(|^) with xi{p) = const for 
isotropic pairing, and xi{p) c< {p^ ~ Py) ^he d-wave 
pairing. 

We have calculated Tc numerically for the 2D super- 
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conductor as a function of magnetic field and spin-orbit 
interaction. The result is shown in Fig. 1 for the magnetic 
field strictly parallel to the surface (to exclude diamag- 
netic currents)!^. 

A 1** order phase transition was initially expected be- 
tween superconducting and normal states, defined by 
comparing their free energies: Fs{T) — Fn — Xn"^, 
which would determine the so-called paramagnetic crit- 
ical field, Hpar {xn - the spin susceptibility in the 
normal phase) . The transition from normal to supercon- 
ducting state is actually (at lower temperatures) a second 
order transition into the LOFF state. The details of the 
phase diagram in the vicinity of the Hpar were studied 
numerically in Rcf.JM. The LOFF phase boundary is 
determined by Eq.(g)or Eq.(U), linearized in Aapiq) at 
an extremal q. The corresponding expression for anoma- 
lous function linear in Aai3{q) is obtained by solving the 
proper Gor'kov equations: 

Fa(3{p, q; iuJn) = -Gi°2(P; iw„)Ap^(p, g)G|j''j (-p+q; -iw„), 

(9) 

where G'^'^^(p; ia;„) is the normal state Green function 

at non-zero hso (p) of Eq. (||) , together with the Zeeman 
term, ^b{c ■ B): 

[iuJn - ^ - hsoip) - IJ-b^B]G'-^^{p; iuJn) = i- (10) 

The spin Hamiltionian on the left side, H — hso{p) + 
fiBO'B, may be easily diagonalized: 



-\~e{p) 
(11) 

(A = ±1 for the two branches, ppy is the y-axis pro- 
jection of the Fermi momentum, and _B || a;) with the 
eigenfunctions, spinors ri^{p) of the form: 



Ae(p) 



(12) 



Substitution of Eq. (^ into Eq. (^) making use of Eq. (||) 
results in a rather cumbersome expression which gener- 
alizes the corresponding Eq.(7) of Ref. jo). We sketch, 
therefore, only a few results for the low-T part of the 
phase diagram in Fig.l. Below we discuss the main 
changes in the shape of the phase diagram, as introduced 
by SO coupling or anisotropy. As for the structure of the 
LOFF phase itself, we assume that the numerical anal- 
ysis done in Ref. ||l4| remains applicable, i.e. the order 
parameter in the LOFF state has the structure of peri- 
odic stripes. 

We state in more detail our results for the limiting 
cases of strong and weak SO interaction, which signifi- 
cantly simplify all calculations. 

a) Strong SO: app ^ ^(0); subsequent analysis, which 
we do not provide here, leads after short calculations to 



our final results: 



q±B, 



\Q\ = 



2fiBH, 



c2 



VF 



tJ.BHc2 = \/2A{0)apF = \l —TcoapF 



2tt 



(13) 



One sees that SO interaction not only enhances the value 
of Hc2 in comparison with the LOFF for 2D model of 
Ref. 1^ , but it also fixes the direction of the structure 
vector. The resulting stripe structure is parallel to the 
magnetic field direction and has the perpendicular space 
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periodicity L = ^i^^^vf- 

According to |T2|, in the limit of strong SO interac- 
tion the spin susceptibility for parallel fields in super- 
conducting state, Xs{T) is non-zero and equal to ^xn- 
This increases the critical paramagnetic field ||l^ only 
by a factor of -\/2, i^sFlpar — ^(0). Comparing this 
with Eq.(p^, one sees that strong SO significantly in- 
creases the area occupied by the LOFF state by a factor 
of \J apFjTca ^ 1. Strong SO scattering by defects also 
enhanses i?c2||l^, but the LOFF state does not exist in 
presence of disorder. Analytical expressions can also be 
obtained at low temperatures. For the dependence of 
the transition temperature on magnetic field, Tc{B), one 
obtains, at T ^ Tcoy/Tco/{apF)- 

TciB) = 5.784 (^^2nTcoapF/7 - ixbB^ l{apFf. (14) 
and 



T,(B) = ^r2,/(27AiBB), T,o^T,ol{apF) « T « T.q 

(15) 

The expression for Tc{B) also simplifies for small mag- 
netic fields near Tcp: 



Teo - 7C(3) (mbS) 



(16) 



Note that while the spin-orbit interaction splits the Fermi 
surface, quasiparticles with the same band spin index 
Eq.(pl]) form Cooper pairs in the superconducting state, 
so that spin-orbit interaction alone does not change Tc- 

b) Weak SO: apF <C A(0); unlike in case of strong SO, 
the Cooper pair is formed mainly by pairing of electrons 
from the FS's with different spin indices. The LOFF 
phase in a 2D superconductor with no SO interaction 
was first analysed in Refs |l^. According to Ref.[|j, 

Msi?c2 = A(0) = V2Hpar, VfQ = 2flBHc2- 

This result is zero order in spin-orbit interaction. The 
direction of q is not fixed with respect to B. Analy- 
sis to the second order in apF results in an anisotropy 
term which again fixes the vector q, as in strong SO case 
Eq.(|l3|), perpendicular to the direction of the magnetic 
field. Indeed, for the critical field as a function of the 
angle /3 between q and B we find: 



HbH,2 = A(0) 



9 9 

a Pf 
2A(0) 



cos /3, 



(17) 
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i.e. the maximum for Hc2 is reached for (3 = ±-. 
For small magnetic fields near we get: 

Too - Te 7C(3) 



47r2 



(18) 



A quadratic dependence on B remains valid for fields 
^lbB ~ apF- However, note the factor two difference be- 
tween Eqs.(p^) and (|l8|): suppression of Tc by magnetic 
field turns out to be slower in case of strong spin-orbit 
interaction than in case when the spin-orbit interaction 
is weak. 

Although a more complicated LOFF periodic super- 
structure is possible, the energy considerations of Ref.^ 
have shown in 3D that the stripe phase is energetically 
more favorable. A detailed study done numerically in 
2D iQ has shown a more complicated than just a si- 
nusoidal shape of the order parameter. We assume that 
the results of Ref. remain valid in all cases considered 
above, so that the LOFF state preserves its striped order 
parameter form. 

The major role of SO is in fixing of the LOFF super- 
structure. Anisotropy fixes the orientation of the LOFF 
stripes as well (see below). A 1** order re-orientation 
transition may be expected corresponding to an abrupt 
change in the direction of the superconducting stripes 
(similar to a spin-flop transition), as the magnetic field 
is rotated in the 2D plane. 

In the above discussion so far we have neglected any 
anisotropy at all. Meanwhile, the anisotropy is, of course, 
important. We address the issue of the effect of pinning 
stripe direction to the particular crystal axis only for the 
d-wave order parameter, since the latter is intrinsically 



anisotropic. Stripes may orient themselves along the di- 
rections of the gap maximums. It can be easily shown 
that Hc2 for the dx2_y2 superconducting order parame- 
ter takes the form: ^bHc2 = ^e^'^^^^co- The critical field 
for the LOFF phase of d-wave is somewhat higher than 
for s-wave (without the SO term Eq.(7)), as supeconduct- 
ing stripes get pinned to the crystal axes by the form of 
the order parameter and the direction of the magnetic 
field. Reorientation transition/twinning is also expected 
for any other cause. 

In summary, we have shown: i) that inhomogeneous 
state in parallel fields extends considerably the low tem- 
perature phase diagram of surface superconductivity 
with increased spin-orbit interaction; ii) all SC charac- 
teristics of the phase diagram in the (i?, T) plane can be 
expressed in terms of Tco, the critical temperature in the 
absence of the field, which is controlled by doping, and, 
hence, not only allows a data comparison in the broad 
range of FET-doping, but may also serve as a method to 
extract the value of the SO interaction; iii) LOFF state 
properties are strongly anisotropic in the plane with re- 
spect to field direction - this, for example, can be seen by 
measuring the anisotropy of the AC susceptibility signal; 
iv) the indispensible feature of the LOFF state must be 
the re-orientation transitions at the field rotation in the 
plane caused by locking of the LOFF order parameter by 
anisotropy. 
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